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Abstract 

Stochastic nonequilibrium exclusion models are treated using a 
real space scaling approach. The method exploits the mapping be- 
tween nonequilibrium and quantum systems, and it is developed to 
accommodate conservation laws and duality symmetries, yielding ex- 
act fixed points for a variety of exclusion models. In addition, it is 
shown how the asymmetric simple exclusion process in one dimension 
can be written in terms of a classical Hamiltonian in two dimensions 
using a Suzuki- Trotter decomposition. 

1 Introduction 

Stochastic models of lattice gas dynamics provide insight into the nonequi- 
librium behaviour of a variety of physical processes, such as surface reactions 
and growth, catalysis and transport phenomena. These models are systems 
of many interacting particles — the dynamics of the particles are prescribed 
in the model definition — and the evolution is typically governed by a Master 
equation. They exhibit steady state phase transitions and very rich dynam- 
ics, but there exists no general framework in which to analyse nonequilibrium 
models. Exact treatments are scarce (see e.g. [HEIIH]) and so approximate 
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techniques are required. To this end we present a scaling treatment designed 
to capture universal and non-universal critical properties of nonequilibrium 
systems. 

The scaling method was developed in detail in 0j. It exploits the well 
known equivalence between the Master equation and the Schrodinger equa- 
tion in imaginary time P - the stochastic lattice gas model is written as 
a quantum spin model. The scaling is achieved using a real space blocking 
procedure [3] in order to thin out the number of degrees of freedom. It has 
been applied to the contact process [HI El, where very accurate results for the 
critical point and certain critical exponents were obtained. Here, we show 
how to adapt the method to models which possess a conservation law (e.g. 
conserved particle number) in order to obtain exact fixed points. Further, 
from a stability analysis of the fixed points, we infer the role of bias in the 
model dynamics. 

2 Quantum Scaling for Exclusion Models 

In the following, we consider exclusion models — models where sites on a lat- 
tice are either occupied by a single particle or vacant. In the quantum formu- 
lation, these models are spin- 1/2 quantum chains. The mapping is achieved 
by interpreting configurations of particles and vacancies in the nonequilib- 
rium model as a configuration of quantum spins, where particles are replaced 
by an up-spin and vacancies are replaced by a down-spin, say. Since the dy- 
namics in the nonequilibrium model become processes involving spin flips, 
they can be expressed in terms of a quantum Hamiltonian. Hence, for a 
lattice containing L sites, the configuration is written |{(7/}) = n^i®! "^ 
where {07} = 01, . . . , o^. We use the notation 07 = +1 to represent an up- 
spin at site I (i.e. a particle in the nonequilibrium system), and 07 = —1 to 
represent a down-spin (i.e. a vacancy in the nonequilibrium system). The 
steady state of the nonequilibrium system is equivalent to the ground state 
of the corresponding quantum problem. 

2.1 Quantum Renormalisation Group Scheme 

We now outline the renormalisation group scheme for quantum systems [3]. 
The first step is to divide the lattice into adjacent blocks, each containing b 
sites, as indicated in figure [TJ The lattice contains V = L/b blocks, labelled 
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Blocking transformation 

v u+ 1 
1 1 

Figure 1: Blocking of the lattice for b = 2. 

v = 1, . . . , L', which will form the sites in the renormalised system. The 
Hamiltonian for the blocked system is written in terms of an intrablock part 
H u , containing all the interactions within block z/, and an interblock part 
H Vil/+ i, containing all the interactions between block v and block v + 1 (we 
assume that there are nearest neighbour interactions only). Hence 

L' 

H = J2Wu + H u , u+1 ] . (1) 

The renormalisation is achieved as follows. We treat the intrablock Hamil- 
tonian H u exactly, and regard the interblock part H UtU+i as a perturbation. 
Thus, we find the eigenvectors of H u and use only those of lowest energy to 
form a truncated basis of states — in a spin-1/2 system we aim to keep the 
two lowest lying eigenstates. The configuration of the renormalised lattice is 
written as a direct product over blocks: 

L' 

iw>=n®i ff ">' ( 2 ) 

v=l 

where {cr u } — <7i, . . . , ay and the block spin states \a v ) are the renormalised 
basis of states obtained from the lowest lying eigenstates of H u . Then the 
renormalised Hamiltonian H' is obtained by writing H' = "Y^ U \H' V + H' uu+1 ]. 
The matrix elements of H' v are given by 

(cr u \H' u \a' u ) = (a u \H u \a' u ) , (3) 

where on the r.h.s., we use the fact that we have assigned eigenvectors of 
H u to the block spin states (since these eigenvectors are orthogonal, i.e. 
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{cr\a') = 5 u>a i although H v is not Hermitian, H' u is diagonal). Similarly, the 
matrix elements of H' u u+l are given by 

( a u\Hl tU+1 \a' u ) = (a v \H^ v+l \a' v ) , (4) 

and thus H' v v+1 contributes the interaction terms in the renormalised Hamil- 
tonian H' . 

Thus we have a prescription whereby we retain only the lowest lying 
eigenstates of a block Hamiltonian in order to thin out the number of de- 
grees of freedom and still retain the features important for criticality — the 
renormalisation is carried out near the ground state of the quantum system, 
equivalent to the steady state of the nonequilibrium system. By assigning this 
truncated basis of states to a block spin variable in an appropriate way, one 
hopes to obtain a rescaled Hamiltonian H' of the same form as the original 
Hamiltonian H but with rescaled parameters. 

2.2 Scaling for the Asymmetric Simple Exclusion Pro- 
cess 

In this section, we apply the above transformation to the asymmetric simple 
exclusion process (ASEP) in one dimension with periodic boundary condi- 
tions. In the ASEP, particles hop to the right (left) with rate p (q) provided 
that the target site is empty. These dynamics can be expressed in terms of 
a quantum Hamiltonian given by 

L 

H = Y,[p ( p i +p i+i - w + i) + ? (prp&i - ^v+i)] • (5) 
i=i 

where = |(1 ± erf) are projection operators and (a^) creates (annihi- 
lates) a particle at site I. Thus the terms crfa^ and a^af +1 generate particle 
hopping to the right and left respectively — they derive from the gain terms 
of the original Master equation (i.e. the terms due to particle hopping which 
increase the probability of finding the system in a particular configuration). 
The terms involving projection operators arise due to the loss terms in the 
Master equation (i.e. the terms which contribute to the probability that the 
system is in a particular configuration provided that no particle performs 
a hop). Therefore probability and particle number (which is related to the 
^-component of the spin) are both conserved. This model also possesses a 
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particle-hole symmetry under interchange of p <-> q. The aim is to maintain 
this duality at all stages of scaling. 

We begin by dividing the lattice into blocks of size 6 = 2. Then we split 
H into an intrablock Hamiltonian H v given by 

Hu = p (p^p-,2 - ^i< 2 ) + q {K,iP„,2 - <i^ 2 ), (6) 

where the suffix u, i indicates that the operator acts at site i in block u, and 
an interblock Hamiltonian H v v+1 given by 

H v , v +\ = P ~ °v,2°Ul,l) + <? (Pu,2 P u+l,l ~ °u,2°v+l,l), (?) 

such that both H v and H„ ;l/+1 possess the duality under interchange of p <-> q 
when + <-> — . 

The next step is to find the lowest lying eigenstates of H u and use them 
to form the renormalised basis of spin states. The ground state of H v is 
three-fold degenerate — since the dynamics conserve particle number H u 
decomposes into b+1 disconnected sectors. This particle conservation can be 
respected in the renormalised Hamiltonian if it is present in the renormalised 
basis of states, therefore the ground eigenstates of H v are organised according 
to their eigenvalue of the block spin operator | Yl\=i a l,i- This leads us to 
define block spin states 

\ + U) = \+u,l,+u,2) , (8) 

|0„) = ^[?l+,,i,-,,2> +pK,i,+,,2>] , (9) 
|-l„) = |-„,i>-„2>, (10) 

for each block v. The corresponding left eigenstates of H v are assigned to 
the left block spin states 

(+1,1 = (+„,!> +„,2| , (11) 
(0,| = (+u,l,~u,2\ + (-,,1,+,,2| , (12) 

(-1,1 = <-„,!, -„ >2 |. (13) 

The left ground eigenstates of a quantum Hamiltonian describing a nonequi- 
librium process are always given by such sums over vectors (where the co- 
efficient of each vector is equal to one) due to conservation of probability. 
By defining the block states in this way, we aim to maintain particle con- 
servation (which we cannot maintain by arranging these states in two linear 
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superpositions forming a spin- 1/2 basis) and the particle-hole duality of the 
model. Thus the renormalised Hamiltonian will describe a spin-1 quantum 
chain. 

We are now able to calculate the matrix elements of H' v and H' uv+ ^. 
Firstly, we note that because our basis states are degenerate eigenstates of 
H v , the contribution due to H' v is a constant and moreover, because the 
groundstate eigenvalue is zero (which is always the case for the ground state 
eigenvalue of quantum systems representing nonequilibrium models) this con- 
stant is zero. The task then is to evaluate the matrix elements of H' vv+1 , 
as prescribed by equation (JJJ). For example, the term in H u<v+ i given by 
a v2 a v+i 1 contributes the matrix elements 

(a u ,a u+1 \H' uu+1 \a' u ,a' u+1 ) = (<T v \a~ 2 \<rl)(a u+1 \<T^ +1>1 \<7 , v+1 ) . (14) 

Thus the operator a~ 2 is replaced by a renormalised operator given by 

\ <— l^Icr~ 2 | + 1„) a~ 2 1 0„) (-l y |o"~ 2 | - 1„) / ^ 

where the 1/ suffix on the matrix represents an operator acting on the trun- 
cated basis of states in block v. Evaluating such matrix elements for all the 
operators appearing in H UtU+ % yields 



<2 


= P +1 + 
v 1 


V pO 
p+q v ' 


P+ - 


p+i 1 q po 

1 v+l T p+q 1 v+li 


P u,2 


_ Q pO 

p+q v 




-^+1,1 = 


P pO 1 p-1 

p+(? - r i/+l "T -^IZ+l) 


<?»,2 = 


^a°' +1 - 

p+q v 


~ U V 1 


a u+l,l — 


p _0,+l , -1,0 
p+ q °v+\ "T* a i/+l > 


°"^2 = 


af'° + - 

V p 


+q U v ' 


°v+\,\ = 


_+l,Q , 9 -0,-1 

a i/+l p+q'Vfl ' 



where P^ projects into the block spin state % and a 1 ^ raises or lowers the spin 
from block spin state i to block spin state j. These expressions are substituted 
into H u>u+ i leading to a renormalised Hamiltonian given by H' = H U)U+ i. 

Thus the renormalised Hamiltonian describes a three-state stochastic pro- 
cess (probability is still conserved) where the z-component of the spin is still 
conserved. In order to obtain scaling equations for the rates p and q this 
Hamiltonian has to be projected onto a basis of spin-1/2 states. To do this, 
we note that raising (lowering) operators in the spin-1/2 system are written 
only in terms of operators which raise (lower) the spin in the spin-1 system. 
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Further, each of the dynamical processes, generated by the raising and low- 
ering operators in the spin-1 system, can be identified with a corresponding 
term involving projection operators due to no transition. Hence we rewrite 

Pv,2 = (1 + ' = (1 + ^tq)Pv+l) 

Pu,2 = (1 + ~^+q)Py 5 ^Wl,l = (1 + £fj)Pv+l7 
<2 = (1 + j&K, ^+1,1 = (1 + ^Kl- 

^ 2 = (i + ^+1,1 = (i + ^y^w 

whence the renormalised Hamiltonian assumes the same form as the original: 
H ' = E b' ( P * +P ^i - + l' (P; p v+i - <tf <Vu)] > (15) 

where the rescaled rates p' and q' are given by 

P' = P( 1 + ^) 2 ' <?' = <?(! + • ( 16 ) 

In order to exploit this renormalisation, we consider the ratio 7 = p/q. 
Stable fixed points for 7 are found at 7* = and 00, and these are separated 
by an unstable fixed point at 7* = 1. Hence symmetric diffusion is unstable 
with respect to bias. At the symmetric fixed point the Hamiltonian (0) 
describes a spin- 1/2 Heisenberg chain whose dynamics are governed by a 
dynamic exponent z = 2. This behaviour therefore is unstable and in the 
presence of any bias the dynamics are described by a new exponent. A Bethe 
ansatz calculation for q = shows that this exponent is z = 3/2 |Hj. That 
we find the exact value for the unstable fixed point is a consequence of our 
preservation of duality at all stages in the blocking. 



2.3 Scaling for the Pair Evaporation and Deposition 
Model 

Another model possessing a duality is a process whereby pairs of particles 
evaporate from adjacent lattice sites with a rate e or are deposed onto adja- 
cent vacancies with a rate 5. The quantum Hamiltonian representing these 
processes takes the form 

H = E i 6 ( P r p m " °i~°ti) + c {PtPfri - ^m)] • (17) 
I 
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The duality in this model is again a particle-hole symmetry under the in- 
terchange of 5 «-> e. There is also a conservation law similar to the particle 
conservation in the ASEP: if we label the sublattice of odd (even) sites by 
A (B), then the density of particles on sublattice A (B) is pa {.Pb)- The 
conserved quantity is the difference in the sublattice densities, pa — Pb- 

The duality of this Hamiltonian is preserved if we divide it into an intra- 
block part given by 

H v = 5 (P- x P- 2 - a+ ff + ) + e (P+P+ - a'^) , (18) 

and an interblock part given by 

H v , v +x = 5 (P~ 2 P v+l,l ~ ^2^+l,l) + 6 ( P ^2 P u+l,l ~ ^2<Vn,l) • (19) 

Again, the ground state of H u is three-fold degenerate with eigenvectors given 
by 

ID = \-v,l,+u,2), (20) 

|2„) = Y^[6\+u,i,+u,2) + e\- v>1 ,- V) 2)} , (21) 
o + e 

|3„) = -„, 2 ) ■ (22) 

We can maintain the conservation of the sublattice densities with respect to 
both the intrablock and the interblock Hamiltonians if we assign spin-1 block 
states in the following way: 

I + If) = 1 1^) j I + lf+i) — , 
|0,) = \2 V ) , \0 V+1 ) = \2 U+1 ) , (23) 

I — lj/) = \3v) , I — lf+i) = |lf+i) • 

Using these states to calculate H' uu+l (H' u is zero as before), and then forcing 
the resulting spin-1 Hamiltonian back into a spin-1/2 basis in the same way 
as was done in the previous section, one obtains a renormalised Hamiltonian 
of the same form as ()17j) with rescaled rates 5' and e' given by 

<5 / = <K 1 + ^) 2 > e' = e(l + ^) 2 . (24) 

Now, in terms of the ratio 7 = 5/e, we again find three fixed points at 7* = 
0, 1 and 00 but now the symmetric fixed point 7* = 1 is stable. Therefore we 
find no dynamic transition in this model — the dynamics are independent 
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of bias. Again, the Hamiltonian (|17|) of the symmetric problem is given by 
that of the spin-1/2 Heisenberg chain. Thus the dynamic exponent is z = 2 
in the biased and unbiased cases. This result is consistent with a mean field 
result which predicts that the diffusion constant is independent of bias, and 
it is supported by numerical simulation 0. 

2.4 Scaling for the Pair Evaporation and Deposition 
Process with Diffusion in One Dimension 

In this section, we apply the renormalisation group transformation to a model 
incorporating both the dynamics of the ASEP and of the pair evaporation 
and deposition process. These dynamics, on a chain with periodic boundary 
conditions, are represented by a quantum Hamiltonian given by 

H = E i 6 ( p i' p iii - a >ti) + e ( p i +p iXi - ^r+i) 

+p {P?Pr + t - W+i) + q {PfPiti - ^V+i)] • (25) 

Exact results have been obtained when 5 = e and p = q in which case the 
model is equivalent to a spin-1/2 XXZ ferromagnet, and also for 5 + e = p + q 
which is the condition that the evolution operator can be written as a free- 
fermion Hamiltonian |2| ■ 

A blocking, with a dilation factor b = 2, is implemented by dividing H 
into an intrablock part H v given by 

+P «i^2 " ^i< 2 ) + Q (Pu,iP^ ~ <i^ 2 ) > (26) 
and an interblock part given by 

H U)U+ i = 5 (P~ 2 Pu+l,l ~ a u,2 a u+l,l) + 6 {P V ,2.Pv+\,l ~ a u,2 a u+l,l) 

+ P (P+2^+1,1 " Wt+l,l) + 1 (P^Pu+1,1 ~ <2^+l,l) • (27) 

The blocking proceeds as in the previous two sections. The ground state of 
H u is two-fold degenerate with eigenvectors denoted 

IU = JTe PI +">i +^ 2 > + e l ~^ --' 2 )] ' (2 8 ) 
|2„> = [q\ +„,i ~u,2) + p\ -„,i +„, 2 >] ■ (29) 
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Effective spin states for the block v are identified by taking the block up-spin 
\+ v ) = \l v ) and the block down-spin \— v ) = \2 U ) for all v. Thus the block 
spin states observe two symmetries of the model: the particle-hole symmetry 
under interchange of rates 5 <-> e and p <-> q, and a symmetry whereby the 
spins on the even sublattice (say) of sites are flipped and the pair evaporation 
and deposition processes are transformed into the hopping processes and vice 
versa. This assignment of block spin states yields a renormalised Hamiltonian 
H' of the same form as but with rescaled rates S',e',p' and q' given by 

£ _ P 3 + <f +pq(5 + e) 
(p + q) 2 
, _ 5e{5 + e + p + q) 
(5 + e) 2 
, _ 5e(p + q) + bp 2 + eq 2 
P ~ (5 + e)(p + q) 

, _ 5e(p + q) + ep 2 + Sq 2 
Q ~ (5 + e)(p + q) 

A flow diagram is obtained where the symmetric fixed point (i.e. where all 
rates are equal) is fully stable. This suggests that the dynamic transition 
in the ASEP is removed when the pair evaporation and deposition processes 
are included — the dynamics in this model are described by the exponent 
z = 2 for all choices of 5, e, p and q (provided S and e are not both equal 
to zero). Also, all the fixed points satisfy the free-fermion condition 5 + e = 
p + q for which exact results are available [3 E] . Moreover, this condition is 
stable (after iterating the transformation an infinite number of times the 
rescaled rates always satisfy the free-fermion condition, regardless of the 
original choice of rates), therefore we do not expect to observe any new 
macroscopic behaviour in the regions of parameter-space to which the exact 
solution does not apply. 



(30) 
(31) 
(32) 
(33) 



3 Trotter Decomposition for the Asymmetric 
Simple Exclusion Process 

One difficulty that arises when constructing renormalisation group trans- 
formations for quantum systems is associated with the non- commutation of 
operators appearing in the quantum Hamiltonian. In this section, we remove 
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this problem by exploiting the Suzuki- Trotter decomposition [TUl HH H21 EE] 
to rewrite the quantum Hamiltonian representing the ASEP as a classical 
Hamiltonian for Ising spin variables in two dimensions. As a byproduct, 
direct contact is made between stochastic nonequilibrium models in one di- 
mension and vertex models in two dimensions [14] . 

The general scheme for the mapping is to split the Hamiltonian into a set 
of operators {Hi} such that 



H = Y / H i , (34) 



i=0 

where each term in Hi commutes with every other (or at least, where any 
non- commutation can be easily dealt with), but [Hi, H m ] ^ for I ^ m. The 
exponential of a sum of operators is then expanded as a product of expo- 
nentials by dealing with the non-commutation through the Trotter formula 



o£i = o H i = li m 



n— >oo 



e n . . . e « 



(35) 



For the ASEP, we begin with the quantum Hamiltonian ( ]5]1 . and divide 
it up in the following way: 

Ho = Ez odd [P p i +p i+i + (1 p r p iXi] > (36) 

Hi = - Ei odd [p *r*S-i + q °t°i+i] , ( 37 ) 
H2 =J2 leyea [p p i +p r + i+^ p r p iXi}^ (as) 

^3 = - Ei even [P ^f^i+l + 9 ^ +<T I+l] ■ (39) 

Ultimately, this choice must be made to reflect the update mechanism - 
here we consider a parallel sublattice update. According to equation ([55]) 
the partition function for the quantum system is written 

Z = lim Z {n) , (40) 

n— »oo 

where we have defined Z( n ) by 



Z (ri) = Tr 



' {3H Q (3H 1 (3H 2 0H 3 



(41) 



The next step is to insert complete sets of basis states into (jUJ). We choose 
the basis \s) = \a±, (72, ... , 0"^), where <Tj is the eigenvalue of the Pauli matrix 
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of; a basis which diagonalises H and H 2 . Hence Zr n ) is expressed as follows 



E 
E 



(si|e' 



.Ml 



s 2 )(s2|e" 



S3) ... (s 2 , 



0^2 /3H"3 

e ™ e n 



nw 



/3Hi 



l S T + l) 



odd 



n 



PH 2 



1 St 6 



l S r+l) 



l«i> 

• (42) 



This insertion of basis states also reflects the update mechanism of the origi- 
nal particle process. Now we interpret the r label as a label for a new spatial 
axis, which we will refer to as the Trotter axis (see Figure EJ). This axis has 
its origin in the time axis of the quantum representation of the Master equa- 
tion. Usually, the Trotter axis is an imaginary time dimension, but since our 
quantum problem is governed by a Schrodinger equation in imaginary time, 
here the Trotter axis represents a real time evolution axis. 

From now, for simplicity, we shall consider the fully ASEP: q = 0. Our 
task now is to evaluate the matrix elements appearing in (|42j) . Since we 
chose our basis states to diagonalise H , the contribution from its exponential 
factor is trivially evaluated. The contribution from the off-diagonal part, Hi 
proceeds as follows: 



E 

.2 odd 



Pp 



n 



] [ (s r |exp 

t odd 

where h[ T is given by 



a i a i+i 



>T+1 



) 



lim exp 

A— >oo 



E h 

Lt odd 



(43) 



In 



(3p 



n 



p+ p- p- p+ 

r l,T r l,T+l r l + l,T r l + l,T+l' 



(44) 



with P± = ±(l±ov) and 07 )T is an Ising spin variable associated with the 
site (l,r), and where 



Lt 



P+ P+ P+ p~ J 

r l,T l,T+l r l + l,T Z + l.T+l ' 
_|_P+ p- p+ p + 

J_P+ p- p- p~ 
* £ l,T £ l,T+i r l+l,T r /+1.T+1 

ID- p- p+ p- 
J_P+ p- p+ p- 



P+ P+ P- P + 

J «,r Lt+1 1 Z+1,t i+l,r+l 

- P~ P + P + P + 

- P~ P + P~ P~ 

- P~ P~ P _ P + 

r l,T r l,T+i r l + l,T r l + l,T+l 

- P~ P + P~ P + 



J-P-P+ P+ P" 

^ r l,T l,T+i r l + l,T l + l.T+V 



(45) 
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The parameter A has been introduced in order to project away unwanted 
configurations, i.e. those which do not represent processes allowed under the 
original stochastic dynamics; for instance, the final term in equation (|4"5)) 
corresponds to a particle hopping to the left, which is forbidden for q = 0. 

To obtain the full effective Hamiltonian for the two-dimensional classical 
system, Zr n ) = limA^oo Tre _//(e ' ff) , we must include the contributions from 
the the Hamiltonians H Q , H 2 and H 3 . Then the effective Hamiltonian can 
be written 

H i,m = £ fcW) + £ ^ (46) 

l,r odd 1,t even 

with 

hif = A$, r - In (?) P;P,; r+1 P^ 1)T Pi 1)T+1 

- In (1 - P;Pi +1 P,; i)T ^; 1)T+1 , (47) 

where the final term here corresponds to the 'stay-put' probability that a 
particle, with a vacancy to its right, does not perform a hop. 

The quantum Hamiltonian (J3J) has now been rewritten as a classical 
Hamiltonian (|4fi)) for an Ising system in two dimensions. The nonequilibrium 
steady state of the original model is characterised by the zero temperature 
behaviour of the quantum model. This limit is recovered in the Ising system 
by taking the extent of the Trotter axis to be infinite. The temperature in 
the quantum system is not the same as the temperature in the classical Ising 
system. Instead, the inverse temperature of the quantum model translates 
into the extent of the classical system in the Trotter direction. Thus the crit- 
ical behaviour in the quantum ground state is expressed through the critical 
behaviour of the finite temperature Ising system in equilibrium. Zero tem- 
perature transitions in the quantum model, which occur as a function of the 
couplings, are caused in the infinite Ising system by varying the temperature. 



3.1 'Plaquette' Structure of Ising Hamiltonian 

The Ising Hamiltonian H^ e ^ representing the ASEP contains a parameter A 
which is infinite. This constraint can be incorporated in a compact fashion 
into a vertex model description. In the representation we have chosen, H^ e ^ 
contains only four-spin interactions. These are shown in Figure 121 as the five 
possible arrangements of spins around a plaquette. 
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UJx U) 2 UJ 3 LOi u 5 

Figure 2: The allowed plaquette configurations with their weights 



T 

T 



Figure 3: Lattice on which the 2d classical Ising model is defined. The 
original spatial axis is labelled by I and r labels the Trotter axis. The config- 
urations around the shaded squares are determined by the allowed plaquette 
configurations shown in Figure 121 



The weights Ui to oo 5 of the allowed vertices are 

1 -f ftp @P 

<jJ\ — UJ 2 — 0J± — 1,^3 = 1 ,^5 = — • (48) 

n n 

All other vertices have zero weight. A lattice configuration is then specified 
by placing the allowed plaquette configurations on the shaded squares of the 
lattice shown in Figure El The open plaquettes in this diagram play the 
role of passive plaquettes, whose configurations are determined only by the 
surrounding configurations shown in Figure 121 Every configuration is allowed 
for an open plaquette and each occurs with weight 1. The partition function 
for this system is 

z= n w «> ( 49 ) 

allowed P^uettes i 

configurations 

As mentioned previously, the division of the quantum Hamiltonian into sums 
over commuting operators must be chosen with a specific update mechanism 
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in mind. From Figure El we see that the division (J36|) to (J39|) combined with 
the insertion of basis states in equation describes a parallel sublattice 
update mechanism — the shaded plaquettes in a row of the diagram represent 
the bonds chosen for an update in one time-step. Alternatively, we may 
have chosen to describe sequential update. This is achieved by dividing 
the quantum Hamiltonian into L pairs of local bond Hamiltonians. One of 
each pair contains only diagonal terms, the other contains only non-diagonal 
terms, and there is one pair for every bond in the system. Then, inserting 
basis states between each pair of local Hamiltonians leads to a sequential 
update mechanism. 

We also note that in the chosen basis, we can write a transfer matrix for 
the evolution of configurations along the Trotter axis. Its matrix elements 
are identical to those appearing in the matrix representation of the discrete 
time Master equation (for a specified update mechanism) with a hopping 
rate given by (3p/n. 

3.2 Scaling 

By rewriting the quantum Hamiltonian as a classical Hamiltonian, the prob- 
lem of non-commutation of operators in the quantum Hamiltonian has been 
removed and the classical Hamiltonian provides Boltzmann weights for phys- 
ical processes (in the sense that each plaquette configuration has a direct 
physical interpretation in terms of the original stochastic dynamics). The 
restricted number of plaquette configurations and the way they are placed 
on the lattice suggest that we can devise direct and straightforward scaling 
procedures. Configurations (and their weights) are built up by piecing pla- 
quettes together in an allowed way. Then one can coarse-grain by matching 
configurations under a change of length scale. This dilation of scale may be 
applied to either the time axis or the spatial axis individually, or to both at 
once. One can also explore how changing the update mechanism effects the 
scaling. 

To illustrate these ideas, consider a putative decimation eliminating sites 
along the time axis for a parallel sublattice update mechanism. This is 
shown, for a particular matching of active and inactive plaquettes, in Figure 
HI The weight for the unrenormalised system is obtained by tracing over 
all allowed plaquettes consistent with a given configuration {o~i}. A scaling 
equation may be obtained, for example, by matching this weight with the 
renormalised weight for a plaquette with spins given by a±, a 2 , cr 4 and a 5 . 
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Figure 4: Decimation of the Trotter lattice along the time axis. Sites labelled 
by a circle are traced over. 

4 Conclusion 

Scaling techniques which are simple to implement have been described within 
the quantum formulation of nonequilibrium exclusion models. In particular, 
the exact fixed points are obtained for the ASEP in which the known depen- 
dence of the dynamics on asymmetry is recovered. For a model involving the 
evaporation and deposition of adjacent pairs exact fixed points are obtained; 
the resulting flow diagram indicates that the dynamics are independent of 
the ratio of evaporation and deposition rates. Further, the stability of the 
free-fermion condition in the model combining the dynamics of the ASEP 
with pair evaporation and deposition indicates that the exact solutions pro- 
vide a complete account of the large scale behaviour. Again, the dynamics 
are found to be independent of bias. 

We note that the projections of spin-1 operators onto a basis of spin- 
1/2 operators in Sections 12.21 and 12.31 are not necessary. Indeed, one could 
continue scaling the system to renormalised bases of higher and higher spin. 
An approach to yield the continuum limit of the ASEP in this way jT^j is 
under investigation. 

We have also shown how to map a quantum Hamiltonian representing a 
nonequilibrium exclusion process onto a classical Hamiltonian in one higher 
dimension. In this way, we have shown how to write the steady state proba- 
bilities for configurations of the nonequilibrium system in terms of classical 
Boltzmann weights. This should enable one to borrow real-space renormal- 
isation group techniques for classical Hamiltonian systems and apply them 
to nonequilibrium systems. 
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